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Abstract. We use spectral decimation to provide formulae for computing the harmonic 
gradients of Laplacian eigenfunctions on the Sierpinski Gasket. These formulae are given 
in terms of special functions that are defined as infinite products. 



1. Introduction 

There are few functions more ubiquitous in Euclidean analysis than the sine, cosine and 
exponential, which are the eigenfunctions of the Laplacian on an interval in R. In the theory 
of analysis on fractals, the Laplacian eigenfunctions arguably have an even more promi- 
nent role, as the Laplacian is the fundamental differential operator on which the analysis 
is based. Despite this, there are a number of interesting open questions about the structure 
of such eigenfunctions. In this paper we consider the local behavior of Laplacian eigen- 
functions on the Sierpinski Gasket {SG), in terms of the harmonic tangents and gradients 
introduced by Teplyaev in [6]. Using the spectral decimation method of Fukashima and 
Shima [ 1 ] (see also Chapter 3 of 0) we give infinite product formulae for the tangents at 
boundary points, and use them to describe the one-sided tangents at junction points. These 
results may be seen as a Sierpinski Gasket version of the well known formulae for the 
derivatives of the sine and cosine functions on an interval, though the precise analogue on 
[0, 1] is more complicated (see Equations l2.2H2.4l ). 

The Sierpinski Gasket is the simplest non-trivial example of a fractal to which the stan- 
dard theory of analysis on fractals applies. We refer to the monographs J2][5l for detailed 
proofs of all results we use from this theory. Recall that SG c R 2 is the attractor of an 
iterated function system consisting of three maps F,-(jc) = (x+qf)/2, where the points qo, q\ 
and q2 are the vertices of an equilateral triangle. This means that SG = UjFj(SG), where 
the sets Fj(SG) are usually referred to as 1 -cells. For a length m word w = w\ . . . w m with 
letters Wj e {0, 1, 2} we define F w = F Wl o ■ ■ ■ o F Wm and call F W (SG) an m-cell. The points 
qi, i = 0, 1,2 are the boundary of SG; the set of boundary points is Vq, and we use V m to 
denote points of the form F w (qi) where w is a word of length m. These V m are vertices of 
the usual graph approximation of S G at scale m, in which vertices x and y are joined by an 
edge (written x ~ m y) if they belong to a common m-cell. Clearly V. = U,„ V,„ is dense in 
SG. 
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The Laplacian A on S G is a renormalized limit of graph Laplacians A m on the m-scale 
graphs: 

Au(x) = - lim 5 m A m u(x) (1.1) 

2 m— >oo 

A m u(x) = ^(«(y) - u(x)) for x e V m \ V (1.2) 

and we say u e dom(A) if the right side of ( 11. Il l converges uniformly on V* \ Vb to a 
continuous function. The function is extended to S G by continuity and density of V*. At a 
boundary point e Vo there is an associated normal derivative defined (with g,+3 = qi) by 

(5 \'« 
- (2u(qd - u{F'l\q i+l )) - u(F?(q i+2 ))). 
„. — i' 

A harmonic function h is one for which Ah = 0, and for any assignment of values on 
Vo there is a unique harmonic function with these boundary values. For this reason we 
identify the harmonic functions with the space of functions on Vo. The harmonic functions 
are also graph harmonic, so it is elementary to compute the values of h on Vi from those 
on Vo, and recursively to obtain the values on V m for any m. It will be useful to formalize 
this by defining the harmonic extension matrices A,, which map the values of h on Vo to 
those on F ; (Vo), by 
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and more generally (h o F w (q ), h o F w (q x ), h o F w (q 2 )) T = A w (h(q ), h(q{), h(q 2 )) T , where 
A w = A Wm A Wm t ■ ■ -A Wl . We usually write this in the compact form 
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The structure of eigenfunctions of the Laplacian is similar to that of the harmonic func- 
tions. Specifically, it is true on SG that if m is sufficiently large then the restriction of 
a function u satisfying —Au — Au from SG to V m gives an eigenfunction of the graph 
Laplacian -A m u = A m u, with 

A m (5 - A m ) = /l m _i (1.3) 
A= 3 -Yim5 m A m . (1.4) 

Note that ( 11.31 ) implies that A m is one of 4(5 + V25 - 4A m -i), but the positive root is only 
permitted to occur for finitely many values of m in order that the limit in ( 11.41 i exists. This 
spectral decimation property was first recognized by Fukashima and Shima (TJ. It is not 
true on all fractals, but on those where it is valid, it gives both a method for computing the 
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spectrum and a recursion for the eigenfunctions [3 1. Let us define 



A (A) = 
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provided A + 2,5. The essence of the spectral decimation method on SG may then be 
summarized in the following theorem, which we have taken from Sections 3.2 and 3.3 
of0. 

Theorem 1.1 (Spectral Decimation Method). If (A + A)u = then there is a sequence 
{Am}m>m a satisfying ( 11.31 ) and ( 1 1 .41 ), and such that A m + 2,5, 6 for m > mo, with the property 
that (A m + A m )u\ v — Oforallm > niQ. The values of u on V m form > mo can be constructed 
recursively using the matrices from d 1 .51 . ( 11.61 ), and ( 11.71 ) as follows. If w = W\W2 ■ ■ -W m 
then 



'■ A Wm (A m )A Km _,(/l m -i) • • • A w (A m(t+1 )u 



Av ) ■ 



(1-8) 



and we call this the spectral decimation relation. 

If A is not a Dirichlet eigenvalue then we may assume mo — 0, at which point the 
condition (A mo + /t mo )w = is taken to be vacuous. The corresponding eigenspace is 

3- dimensional and parametrized by the values of u on Vo. 

If A is Dirichlet several possibilities occur. We indicate the initial configurations, all of 
which may then be continued by the spectral decimation formula. A spanning set for the 
configurations when mo — 1 are shown in Figure\I] The one on the left has Ai — 2 while 
those on the right have A\ — 5. If mo > 2 then A mo — 5 or A mo — 6, and in the latter case 
A/mo+i — 3. Those with A mo — 5 are formed from scaled and rotated copies of the functions 
on the right in Figure\l\ arranged so that their normal derivatives cancel. A basis of chains 
for mo — 2 is shown in Figure^ those for general m are naturally indexed by the loops in 
V m , plus two strands connecting points of Vo- In the case A mo — 6 the eigenfunctions are 
indexed by points in V mo -i \ Vo; a basis is obtained by scaling and rotating two copies of 
the function on the left in Figure\3\and gluing them at the chosen point, as shown on the 
right in Figure\3]for the case mo — 2 and a point in V\. 
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2 -1 0-12-10 

Figure 3. Eigenfunction construction in the case A„, = 6. 

2. Tangents to Eigenfunctions 

In (6|, Teplyaev introduced the notion of a harmonic gradient and harmonic tangent for 
functions on the Sierpinski Gasket. He also proved that functions in the domain of the 
Laplacian have harmonic gradients at all points in a set of full measure, and that functions 
with Holder continuous Laplacian have harmonic gradients at all points. Eigenfunctions 
of the Laplacian fall into the latter category because of the well-known fact that continuity 
of A 2 u implies Holder continuity of Au (El, Lemma 2.2.5). 

Deflntion2.1 ([6|, Section 3). Let w — W\ W2 ■ ■ ■ be an infinite word, and [w] m = wyWz ■ . ■ w„ 
be the length m truncation of w. If u is a function on S G we let H[ W ] m u be the harmonic 
function on SG which coincides with u on F[ W ] m (Vo), so 

H [w]m u = u\ F[w]in {Vo) = K\ ■ ■ ■ A wi u |f Mm ( v ) ■ ( 2 - 1 ) 

Define the harmonic tangent T w u of u at w to be lim^oo H[ W ] m u, if the limit exists. It should 
be remarked that there can be two words w and w' such that F W (SG) = F W >{SG). We will 
nonetheless treat the tangents T w and T K > separately, as they are rarely equal. 

The harmonic gradient is defined in the same way, but using the space of harmonic 
functions with average zero and the projection of the action of the matrices A, to this 
subspace. It is evident that if u is continuous then the gradient exists whenever the tangent 
exists, and conversely. See [4J for details. 

If we consider the interval [0, 1] rather than the Sierpinski Gasket then it is clear that the 
harmonic tangent at jco is the vector (L(0), L{\)) characterizing the unique linear function 
L(x) having the properties L(xo) = /(xq) and L'{xq) = f'(xo). For an eigenvalue A e R 
which is not equal to ji 2 k 2 for any k e N, it may readily be verified that the harmonic 
tangent is given by 

//(0)\ 



M(x )^\ (2.2) 
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where 

M(xo ) = ^-l a + b h a ~ b \ (2.3) 
sin VI I" " b a + b l 

and 

a = sin ((1 - x ) VI) , b = VAx cos ((1 - x Q ) VI) . (2.4) 

We cannot obtain as explicit a description of the harmonic tangent on SG; however, we 
produce formulae that permit its computation at any point of V». The key observation is 
that when u is an eigenfunction, the computation of the gradient has a particularly elegant 
structure. Recall from Theorem 1 1.1 1 that the values of u on f[ W ] m (Vo) ma y be computed 
using the spectral decimation method, meaning that starting from a scale niQ they can be 
obtained as in (11.81 l. Combining this with ( 12.11 ) we see that 

T w u= lim A~]---A~] ■ A Wm (A m ) ■ ■ ■ A w (A m+1 )u\ lv , 

= « ■ ' <Jfc<+, ■ -< ■ A ^ ■ ■■A w ,„ 0+ M mo+1 )) u\ Pwj ^ (2.5) 

in which we know the limit exists by Theorem 3 of [ 6 1 . A special case occurs when F W (S G) 
is a point in V*, because in this case, all but finitely many letters in the word w are equal to 
a single letter i. By taking mo to be sufficiently large we see that it is useful to understand 
the limit 

lim A~ k ■ Ai{A m+k ) ■ ■ ■ Ai(A m+l ) 

and it is evident from the symmetry of the matrices A, and A,(/l) that it suffices to deal with 
the case i = 0. 

Theorem 2.2. Let a = (0, 1, 1) T , /3 = (0, 1, -if, y m = (4,4 - A m ,4 - A m ) T . If neither of 
the values 2 or 5 occur in the sequence {A m } m>m , then 



4/1 °° A 

lim A "* • A (A mo+k ) ■ ■ ■A (A mo+1 )a= - ■ _ — Q(l - Jfl)a 



, 2A 
lim A * • Aa(A mo+k ) ■ • ■ A (i mo+ i)yS = /? 

>°o 3 ■ j /t mo 

lim Aq* • AoOJ^+t) ■ • -AoUmo+Ormo = (4,4,4) r 

Proof. Observe that a and/? are eigenvectors of Ao(A), with eigenvalues (6-/l)(2-/l) _1 (5 — 
A)~ l and (5 - /I) -1 respectively, provided A + 2,5, As Aq - Ao(0) is a special case, we 
compute immediately that 



lim A k ■ A (A„, 0+k ) ■ ■ ■ A {A mo+l )l3 = lim 5* f~[(5 ~ 'W;) l P- 

However induction on ( 11.3l l implies that = /i mo+/ t nj=i(5 - Am) and therefore 

2A 

lim A * ■ A Q (A mo+k ) ■ ■■AoiA^+OP = lim 5 k A n , 0+k A-Jj3 = 
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The corresponding computation for a can be simplified by observing from (11.3b that (6 
A m ) = (3 - A„, + i)(2 - A m +\), whereupon 



lim A * • A (A ma+k ) ■ ■ ■ A () (A mo+l )a 



iW£ff}_ 
V 3 ' 1 1 rs — 



(6 - /l mo + ;) 



y =1 (5 — A„, B+ j)(2 — A mo+ j) 



%-il^)(^nd-)( 3 - t n(3-w 1 )). 

= f ? V ^ )rt(i - 

\2-A mo+1 )\3-5">°A mo nj\ 3 > 



where we used the previously computed limit for the middle factor, and the fact that A„ 
0(5~ m ) as m -> oo from ( fTH i. 

For y m the situation is a little different. Observe that 



Ao(A m+ \)y„, - (2 - A m+ i) '(5 - A m+ \) 1 



4(2 - A m+ i)(5 - i m+1 ) 
4(4 - /l m+ i) + (4 - ,i m )(6 - A m+l ) 
4(4 - /l m+ i) + (4 - A m )(6 - A m+i \ 

however we can perform the following simplification from ( 11.31 : 

4(4 - A m+l ) + (4 - A m )(6 - A m+ i) = 4(4 - A m+1 ) + (4 - 5A m+l + A 2 m+l )(6 - A m+X ) 

= 4(4 - A m+{ ) + (4 - A m+ i)(l - A m+ i)(6 - A m+l ) 

= (4-A m+ i)(10-7^ m+ i+^ +1 ) 
= (4 - /l m+ i)(5 - /l m+ i)(2 - A m+ \). 

Inserting this into the previous computation shows that Ao(A m+ i)y m = y m+ \ provided A m 
2,5, and therefore A (A mo+k ) ■ ■ ■ A (A m(j+ i)y mo = y mo+k . To proceed we must apply A~ k to 
Tmo+k, which is most easily done by writing it in terms of eigenvectors as y,„ a +k = (4, 4, 4) T - 
A mo+k a. The result is 



]imA Q k ■A (A mo + k )---A (A m<l+1 )y mo = lim A () *((4,4,4) r - A ma+k a) 

k— >oo k— >oo * ' 

= lim((4,4,4) r -(^)^ mo+i ) 



(4,4,4) r 



because A m = 0(5 m ) as m — > oo. 



Theorem 2.3. Suppose w is a word of the form w — [w] ko 000 ■ ■ ■ , that (A + /1)m = 0, and 
that niQ is chosen large enough that the spectral decimation formula holds with A m £2,5,6 
when m > otq. Ifk — maxj&o, mo] then 



T K u = (A'] ■ ■■A-l)M (A,k)A Wk (A k ) ■ ■ ■ A Wmo+1 (A m+1 ) 



n(Vb) 



where 



M (A,k) 



1 

(4-A k )A Tk (A) A(2 Tk (A)+l) A(2 Tk (A)-l) 

3-5*A 3-5*A 3-5*4 

I _ (4-A k )Ar k (A) A(2T k (A)-l) A(2T k (A)+l) 



1 - 



3-5*4 



3-5 k A k 



3-5*4 



(2.6) 



(2.7) 
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and 

J— - 



,,..,, t (V ) 



Proof. From ( 12.51 ) we have 

7> = (A-[ ■ ■ • A^)(Hm V"~ k) ' A o0« ' ' ' AoMn-i)) « 

= ■ • OjMm^ ' A„(i m ) ■ ■ • A () (^ + i)) A lv ,(4) ■ ■ ■ A H . mo+1 (^ mo+1 ) ^ 

because w ; - = for j > k and spectral decimation applies for j > mo- The result is therefore 
equivalent to 

Mo = lim A (m - k) ■ A (A m ) ■ ■■A (A k+1 ) 
which follows from Theorem l2.2l □ 

Remark 2.4. The function Tk(A) may appear to depend on the sequence {Aj}, but in fact this 
sequence is uniquely determined by A. Indeed, there is an entire analytic function *P(z) with 
the property that Aj = T(5"U). To see this, let ij/{z) = z(5-z) and ifrJz) = i// om (p' m z). The 
sequence if/ m (z) consists of entire functions with iff„,(0) = and i//' m (0) = I, so is normal 
with limit ^(z), a power series for which may be computed recursively. It is then clear that 
( T(5 _ -'/i) = lim„, *P(|5 m_7 /l m ) = lim,„ if/ m -j(A m ) = Aj. Moreover, we may define 

r ( A) = L^fT (l _Z^) (2.8, 

at which point t^(/1) = T(5"*/l). In the same way that there are special functions associated 
to differential equations in Euclidean analysis, we suggest that the functions ^(A) and Y(/l) 
should be considered to be special functions in analysis on the Sierpinski Gasket. In terms 
of these functions, ( 12.7b has the form 

1 

i(4->F(5- 1: i))r(5- 4 ,l) A(2T(5- l A)+l) A(2T(5- k A)-l) 



M (A,k) = 



1 3-5 k W(5- t A) i-5<"¥(5- k A) 3-5 tl P(5-*/l) 

i(4-"P(5^i))T(5-*i) /)(2T(5-*/i)-l) ,l(2T(5- t /l)+l) 



As a particular consequence we may compute the normal derivatives of the eigenfunc- 
tions at points of Vo, because they are the same as the normal derivatives of the tangent 
functions. We expect this observation to have applications in the construction of a resol- 
vent for the Laplacian. 

Corollary 2.5. If (A + A)u — and the spectral decimation formula holds with A m + 2, 5, 6 
for m > 0, then the normal derivative of u at qo is 

d„u(q ) = ((4 - Ao)u(q ) - 2u{q{) - 2u(q 2 )) — ■ 

Excluding the values 2,5 and 6 from the sequence \A m } in Theorems 12.21 and 12.31 is 
necessary because they occur precisely in the Dirichlet case, where the boundary data 
vanishes and cannot be used to determine the tangent. Nonetheless, Theorem 12 . 3 1 may be 
applied to find tangents to Dirichlet eigenfunctions in a simple fashion. The reason is that 
the description of the Dirichlet eigenfunctions given in Theorem 11.11 ensures that we need 
only compute the harmonic tangents of the functions in Figure [TJ and the left of Figure [3] 
All harmonic tangents to Dirichlet eigenfunctions are then obtained from these by scaling 
and taking suitable linear combinations. 
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For the basic element used to construct the 6-series (Figure we can directly apply 
Theorem 12 . 3 1 with Ai = 6 and Ai = 3. For example, if the top vertex in Figure[3]is qo and 
w = • • • we have 

(0^ 



A 

T w u = - 



1 

-1J 



To calculate the harmonic tangents of the basic element of the 2-series (on the left in 
Figure [TJ we apply Theorem 12.31 to the function shown at left in Figure [4] starting the 
spectral decimation at each of the values A\ = 5± f^ . For the 5-series there are two basic 
elements (shown at right in Figure [TJ and we proceed by calculating harmonic tangents 
for the initial configurations shown in the center and on the right of Figure @] starting with 
A] - ^fi. The harmonic tangents of all 2 and 5-series eigenfunctions then coincide with 
scaled and rotated copies of these pieces and their negatives, assembled in the obvious 
manner. 
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